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Abstract 

We establish nearly optimal rates of convergence to self-similar solutions of Smolu- 
chowski's coagulation equation with kernels K = 2, x + y, and xy. The method is 
a simple analogue of the Berry-Esseen theorem in classical probability and requires 
minimal assumptions on the initial data, namely that of an extra finite moment con- 
dition. For each kernel it is shown that the convergence rate is achieved in the case of 
monodispcrsc initial data. 
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1 Introduction 

Smoluchowski's coagulation equation 



2 



d t n(t, x) = - [ K(x - y, y)n(t, x - y)n(t, y)dy - I K(x, y)n(t, x)n(t, y)dy 
lo Jo 



is a fundamental mean-field model for cluster growth that arises in a wide range of fields, 
including physical chemistry, astrophysics, and the dynamics of biological systems (see [2] 
for a review). Here, n(t,x) is the density of the number distribution of clusters of size 
x E (0, oo) at time t > and K(x, y) is a symmetric rate kernel. In the case of discrete sizes 
/ G {1, 2, • • • } an analogous equation for the coefficients n/(t) of the number distribution is 



1 

(1.2) d t ni(t) = 2 XI K i-jJ n i-j( t ) n j( t ) ~ X K i,i n «( t ) n i( t )> 

3=1 3=1 

with Kij = K(l,j). The continuous and discrete cases can be considered together via 
the weak formulation of Smoluchowski's coagulation equation, given in terms of a moment 
identity for the number distribution n(t,dx): 

(1.3) d t I <f>(x)n(t,dx) = - [ [ (<p(x+y)-<j>(y)-<f>(x))K(x,y)n(t,dy)n(t,dx), 



(0,oo) 2 J(o,oo) J(0,oo) 
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where <j> is a suitable test function. We direct the interested reader to [6, 7] for further dis- 
cussion on the well-posedness and asymptotic behavior of measure- valued solutions n(t, dx) 
to (1.3). A general survey of existing literature in stochastic coalescence is given in [2], with 
more recent work reviewed in [8] . Note that throughout this paper we try to keep the same 
notation as in [7]. 

The present work is restricted to the homogeneous, 'solvable' kernels K = 2, x + y, xy, 
for which K(ax, ay) = a Y K(x, y) with 7 = 0, 1, 2. It has been shown in [6] that (1.3) admits 
a one-parameter family of self-similar solutions whose domains of attraction under dynamic 
scaling are characterized by the tails of the initial data. In particular, there is a unique 
self-similar solution with finite (7 + l)th moment, which has exponentially decaying tails 
and attracts all initial data that satisfy this finite moment condition. In [7], it was shown 
that with an additional integrability hypothesis on the Fourier transform of the initial data 
one has uniform convergence of densities (in the continuous case) or coefficients (in the 
discrete case) to a self-similar solution as t approaches the time horizon T 7 , where T 7 = 00 
for 7 = 0, 1 and T 7 = 1 for 7 = 2. The proof of this result is similar to that of uniform 
convergence of densities in the central limit theorem (see Feller [5], Section XV. 5). 

Here we establish near-optimal L°°-rates of convergence to the exponentially decaying 
self-similar solutions of Smoluchowski's coagulation equation with K = 2, x + y, xy. Con- 
tinuing the analogy with the classical central limit theorm (CLT), this result corresponds 
to the Berry-Esseen theorem for rates of convergence to the normal law ([5], XVI. 5). The 
method is simple and robust: it holds for all initial distributions, requiring only an addi- 
tional finite moment condition beyond those needed for well-posedness and convergence to 
the self-similar profile. In particular, it is true when the initial distribution has a density 
or is a lattice measure. Our work improves upon recent results of Canizo, Mischler, and 
Mouhot [3] as it holds for a large class of initial data and for all of the solvable kernels. 
More broadly, it once again demonstrates the utility of applying methods from classical 
probability to study asymptotic behavior under rescaling for integral equations of convolu- 
tion type. Additional estimates that make use of the present work, such as large deviations 
theorems, will be developed elsewhere. 

We present the results for the solvable kernels in a unified framework. With to = 1 for 
7 = 0, to = for 7 = 1, 2, and no(dx) = n(to, dx), define the moments fij = f, Q ^ x^no(dx) 
of the initial data. Assume that the 7th and (7 + l)th moments are finite, for which we 
scale x and no so that /U 7 = /i 7 +i = 1. As shown in [7], finiteness of the 7th moment ensures 
well-posedness of (1.3), while finiteness of the (7 + l)th moment guarantees that the initial 
data are in the domain of attraction of the self-similar solution with exponential decaying 
tails. These self-similar solutions are absolutely continuous and are explicitly given in terms 
of their densities (see [2, 6, 7]) by 



(1.4) 




Here, the profiles n* i7 for x > take the form 





—X 



xh^i(x) = x 2 n*^(x) 




1 



x/2 



= e 
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and the time-dependent moments m 7 and scalings A 7 are 

(1.6) m (t) = t-\ m 1 {t) = \, m 2 (t) = (1 - t)" 1 



(1.7) Xo(t) = t, X 1 (t) = e M , \ 2 (t) = (l-t)- 2 . 

The main result of our work is as follows. Define the time parameter r 7 (i) = f£ rn-y(s)ds, 
explicitly given by 

(1.8) r (t)=logt, n(t)=t, T 2 {t) = \og{l-ty 1 . 

If we assume finiteness of the (7 + 2)th moment of the initial data, we have exponentially 
fast convergence to self-similar form: 

Theorem 1.1. Let uq be a positive measure with Jj Q ^ x' y no(dx) = Jj Q ^ x^ +1 riQ{dx) = 1. 
Suppose that the additional finite moment assumption fr 0oo \ x 7+2 no(dx) < 00 holds and let 
n(t,dx) be the measure-valued solution to Smoluchow ski's equation with K = 2, x + y, or 
xy and initial data no(dx). With rescaled solution 

A (i) 7 

(1.9) n(r 7 ,d£) := 7 n(t, Xy(t)dx) 

777"-y ^ V J 

define the corresponding distribution functions 



(1.10) F 7 (r 7 ,x)= / y 7 n(r 7 ,(iy), F^(x) = / y 7 n* i7 (y)<iy. 

J(Q,x] J(o,x] 

Then for r 7 G [0, 00), 

sup|F 7 (T 7 ,i) - F* i7 (x)| < C(/i 7+2 )(l + r 7 )e~ T7 

where C(// 7 +2) is a constant that depends only on /i 7 +2- 

Since we are primarily concerned with an asymptotic rate of decay to self-similarity, 
we make no effort to optimize the constant in front of the exponential term. We also 
demonstrate that exponential decay rates are achieved in the case of monodisperse initial 
data no(dx) = 5\{dx). Specifically, we show that for these initial conditions 

sup|F 7 (r 7 ,x) - F* )7 (z)| = <3(e" T7 ). 

The method of proof for these results requires us to work in the Fourier domain. To do 
so we make use of a smoothing argument given by Feller [5], XVI. 3, which for completeness 
we restate here without proof. The distribution functions in (1.10) satisfy the assumptions 
of the lemma. 
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Lemma 1.2. Let F and be probability distribution functions and assume \F^(x)\ < 1. 
Consider the Fourier transforms of their corresponding measures, 



u(ik) 



f e~ ikx F{dx), u*(ik)= f e~ ikx F*(dx). 
J (—00,00) J (—00,00) 



(—00,00) 



Assume that u'(0) = it*(0) — that is, equality of the first moment of F(dx) and F*(dx). 
With the mollifier 

1 — cos(Tx) 



tp T {x) 



ttTx 2 



define 



A(x) = F(x)-F,(x), A T (x) = 2/j T (x) * A. 



Then 
(1.11) 

X 

In terms of Fourier transforms, (1.11) is 



24 

sup|A(x)| < 2sup|A T (x)| + — . 

x x "K-L 



(1.12) sup \F(x) - F*(x)| < - sup 

x 7T x 



T ^ikx 



ik 



(u(ik) — u*(ik))dk 



+ 



24 

7rT' 



In the following sections, we will need to approximate Fourier transforms near the origin 
by their Taylor expansions. To do so, we will use the following basic result also due to 
Feller (see [5], XV. 4 and the subsequent appendix). Let C + = {z £ C : Re(z) > 0} and 
C_ = {z G C : Re(z) < 0} be the open left- and right-half of the complex plane, and denote 
by C+ and C_ their respective closures. Let F{dx) be a probability measure supported on 
[0,oo), which is therefore characterized by its Laplace transform. By analytic continuation 
we can define its Fourier-Laplace transform 



;i.i3) 



it s 



[0,00) 



e- sx F{dx), s G 



which is approximated as follows: 

Lemma 1.3. Suppose F(dx) is a probability measure on [0, 00) with finite 'yth moment 



(1.14) 



M 7 := 



[0,00) 



x^Fidx) < 00, 



where 7 > 1 is an integer. Then for £, s G C+, the Fourier- Laplace transform (1.13) satisfies 



(1.15) 



u(£ + s)-\ u(0 



1! 



-* + ■ ■■ + (-!)' 



( 7 -l)! 



<^i-r- 
7! 



Proof. We begin by showing that for s G C_, 
(1.16) 



s s 7-1 



< 



7! 
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Define gi(s) = Jq e z dz and g 7 (s) = 7 7 _i(z)cfa, where the integrals are taken over the 
straight line segment beginning at the origin and ending and s. A simple calculation shows 
that g 1 (s) = e s - (l + j, + • • • + 
note that 



7— i \ 

( 7 -i)! J • Then (1.16) follows by induction. To see this, 
•1 



\gi(*)\ 

since |e rs | < 1. Assuming |<7 7 _i(s)| < ( 7-1 ) 



s / e TS dr 


17-1 



< s e 



we obtain the desired estimate: 



b 7 ( s )l 



i / 5 7 _i(rs)dr 




< 



(7-l)!7o 



7! 



Equation (1.15) now follows directly. Using (1.16) and the definition (1.14) yields 



it(f + s) 



—s + 



1! 



+ (-1)" 



,7-1 



[0,00) 



sx 

IT 



+ ••• + 



i) 



( 7 -l)! 

'sx) 7-1 



7 -l 



(7-1)! 



<r M M 17 
< — f\s\ '. 

7! 



□ 



Taking £ = in (1.15) and using that M 7 
expression 



tt^(O), gives the particularly useful 



(1.17) 



1 



<^i.r, 

7! 



s 6 



Let us outline the remainder of the paper. We first prove Theorem 1.1 in the case of 
constant coagulation kernel in Section 2, where the method of proof is simpler and more 
easily applied. The same method is then used in Section 3 to establish the result for the 
additive kernel. Finally, the result for the multiplicative kernel is given in Section 4 by 
transforming to the case of additive kernel through a well-known change of variables. 



2 Rate of convergence for the constant kernel K = 2 

In this section we restate and prove Theorem 1.1 for the constant kernel case. That is, we 
show: 

Theorem 2.1. Let uq be a positive measure such that f, Q ^ rio(dx) = f, sXno(dx) = 1, 
and fj,2 = /(Qoo) x2n o{dx) < oo. With r = logi and in terms of the rescaling (1.6-1.7), let 
n{r,dx) be the rescaled solution (1.9) to Smoluchowski's equation with K = 2 and initial 
data no(dx). Then with 



and t G [0, oo) , 
(2.1) 



F(t,x)= / n(T,dy), F*(x) = / h* fi (y)dy. 

J(0,x] J{0,x] 



sup|F(r,f) -F t (x)\ < C(/i 2 )(l + r)e- 
£>0 



where C(fi2) is a constant that depends only on fj.2- 
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2.1 Proof for the constant kernel 



The proof is can be summarized as follows. For s £ C+ define the Fourier-Laplace transforms 
of the rescaled number measures as 



u(t,s)= / e sx h(T,dx), u*(s) = / e sx n*(x)dx 
J {0,00) J {0,00) 

and let uq{s) = u(0,s). By (1.12), we have that 

24 



(2.2) sup|F(r,x) - F*(x)\ < -sup 

x>0 x>0 



(u(r, <r) — u*(a))da 

■iT a 



The uniform rate of convergence of the distribution function is thus given by the convergence 
rate of u(t, s) to tt*(s). As shown in [6, 7], u satisfies a first-order PDE that can be solved 
in terms of u$ by the method of characteristics in C+. The finiteness condition on ^2 is a 
statement about the decay of the tail of no, and hence, about the regularity of uq near the 
origin. Using (1.17) to estimate the difference uq(s) — tt*(s) near s = 0, we can propagate 
this estimate outward at the growth rate of the characteristics. Combining this with Lemma 
1.2 gives us the desired result. 

1. Evolution of characteristics: As shown in [6], the Fourier-Laplace transform u(t,s) 
solves the PDE 

(2.3) d T u + sd s u = — u(l — u). 

This equation can be derived directly from (1.3), or from (1.1) if the initial number measure 
has a density as in [7]. Fix an initial point so € C+. The solution to (2.3) is given in terms 
of characteristics s(t;0, sq) starting from sq at r = 0. The characteristics satisfy 

ds 

(2.4) — = s, s(0; 0, so) = so £ C+, 
whose solution is 

(2.5) s(t;0,s ) = e T s G C+. 

Note that the characteristics are independent of the initial data uq. Geometrically, they 
are rays emanating from the origin at speed e T . In particular, the map so > s leaves the 
imaginary axis invariant. Along characteristics we have 

du 

— = - u (l - u), 

which can be integrated to yield the solution in terms of so: 

uq(s ) 



(2.6) u(r,s) 



l-Uo(«o)(l-e- T ) 
The Fourier-Laplace transform of the self-similar solution, 

(2.7) u,(s) = 



6 



can also be expressed in terms of characteristics. Fixing Sq £ C+, the characteristics 
s*(t;0, Sq) corresponding to the solution of (2.3) with initial data u* are 

(2.8) s*(r;0,s* ) = e T s* eC + . 

These characteristics are identical to (2.5), as expected. Since u*(s) is a time-independent 
solution to (2.3), 

(2.9) «,(s)=e- T Ms * o) 



l-n,(s*)(l-e-^)' 

2. Estimates near the origin: Now we utilize the additional finite moment assumption 



fJ-2 

and that 



x ho(dx) = / x 2 no(dx) < oo 

(0,oo) J(0,oo) 



2 = / x 2 n*{x)dx = I x 2 n*(x)dx 
«/(0,oo) J(fl,oo) 

to obtain approximations for u$ and it*, respectively, near the origin. By (1.17), 

(2.10) \u Q {s) - {1 - s)\ < ^|s| 2 , M«)-(l-«)| < |s| 2 , sGC+. 

We can also estimate u(r, s) near the origin using the time-dependent moment 
7712(7") = / x 2 h(r,dx) = / x 2 n(e T ,dx). 

J(0,oo) V(0,oo) 

From (1.3) (or (1.1) if n(t,dx) has a density) it follows that $tm 2 (r) = 2 — ^(r), so 

m 2 (r) = 2 + (/i 2 -2)e^ T . 
Therefore, (1.17) implies that for all r > 0, 

(2.11) Kr,s)-(l- S )|<^M 2 <f | S | 2 , S GC + . 

3. Propagation of estimates: Consider the r.h.s. of (2.2) with T = 5e T , where 5 = 
3 (^2) > is a length scale characterizing the region about the origin where the difference 
between uq and it* is small by (2.10). Explicitly, define 



(2.12) 2S(f, 2 ) = ^1 + 2(1 + 

so that (1 -I- 1^2/2) 5 = (1 + S)~ 1 /2. As expected, if /i 2 — > 00 then 5-^0. 

4- Backward characteristics: Denote the common set of characteristics for u and it* by 
cr(r; 0, <7o) with o"o G C+. Now define the backward characteristic ctq{t, ao) as the inverse of 
the forward mapping ao >->■ <t(t; 0, ao). That is, for 10 € C+, 

(2.13) o~o(t, uj) = c(0; r, cj) 
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Also note that by (2.5) or (2.8), the Jacobian of the transformation <jq *— > a is 

da 



(2.14) 



da 



We now estimate (2.2) by considering separately the integral near and away from the origin. 
4- Bounds for < \o~\ < 5e~ T : Applying (2.10) and (2.11), it is simple to see that 



/ (u(r,a) — u*(a))da 

Jo<|cr|<<5e- T a 



(2.15) 



i — :\u(r,a) — u*(a)\d\a\ 
'o l°1 

<(l + f)^. 



5. Bounds for be T < \a\ < 5e T : Substituting (2.6), (2.9) into the r.h.s. of (2.2) and 
changing variables to ctq with (2.5) and (2.14) we have that 



(2.16) 



(u(t, a) — u*(a))da 

8e- T <\cr\<8e T a 

„ct(t;0,(To)x 



{u (a ) - U*(<7o)) 



e-^<\a \<8 °"( T ; 0,cr ) (1 - u«(0o)(l - e~ r ))(l - U ((J )(1 - e- T )) 



d(T 



Since ctq is pure imaginary in the region of integration, (2.7) implies 



|1 - u*(a )(l - e T ) 



e T + a 



l + o-Q 



> 



l + <5 



and by (2.10), 



1 kol 



|1 - uo(ao)(l - Ol > |1 - u*(a )(l - e" T )| - (l + y) K| 2 > ^ ■ { 



+ 5 



Therefore, 



/ 

J 5, 



(u {a ) - u*(cr )) 



Se-^<\a Q \<8 ct(t; O.o-q) (1 - u* (<t ) (1 - e T ))(l - u (<7o)(l - e r )) 



< 4 (l + ^)(l + 5) 2 e- ^d\a \ 



(2.17) 



8(l + f) (l + 5) 2 re- 



5. -Rate of convergence: Substituting (2.15) and (2.17) into (2.2), we obtain the desired 
result (2.1). 



2.2 Near-optimality of the rate for K = 2 

We demonstrate that the exponential rate is achieved for monodisperse initial data. Using 
the explicit solution given in [2] the solution with K = 2 and initial data no(dx) is 



(2.18) 



n(t,dx) = J]V 2 ( 1 - - 5j(dx), t€[l,oo). 
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Here, 6j(dx) is the Dirac measure centered at j £ N. In similarity variables (1.9) the 
distribution function for n(r, x) corresponding to (2.18) is 

F(r,x) = e- T Yl (l-e- T ) j ~ 1 = l-(l-e" T ) rer:J;1 " 1 

l<j<e T x 

where \x] = min{j € N : j > x}. The limiting distribution is 

F*(x) = 1-e - *. 

We have that 

(2.19) F(t, x) - F*(x) = e~ x (l - e a ^ T ' x ^ = e~ x (a(r, x) + 0(a 2 (r, x))) 
where 

(2.20) a(r, x) = x + ( ^x] - 1) log(l - e~ T ) = o((l + ^\ e~ T ^j 

uniformly in x as r — > oo. Since < (1 + x/2)e~ x < 1, substituting (2.20) in (2.19) gives 
the desired result: 

sup|F(r,x)-F*(x)| =0(0- 

x>0 



3 Rate of convergence for the additive kernel K = x + y 

We now consider the case of the additive kernel: 

Theorem 3.1. Let uq be a positive measure, Jj Q ^ xn$(dx) = Jj Q ^ x 2 no(dx) = 1, and 
Us = f^ Q ^ x 3 no(dx) < oo. With h(t,dx) the rescaled solution to Smoluchowski's equation 
with K = x + y and initial data no(dx) let 

F(t,x)= / yn(t,dy), F*(x) = / yn*,i(y)dy. 

J(0,x] J(0,x] 

Then for t 6 [0, oo), 

(3.1) sup|F(i,x)-F*(x)| <C(/z 3 )(l + t)e _ * 

Z>0 

where C(fis) is a constant that depends only on ^3. 



3.1 Proof for the additive kernel 

The analysis for the additive kernel is identical in spirit to that of the constant kernel but 
is complicated by the fact that the characteristics of the Fourier-Laplace transform are no 
longer rays. We use a contour deformation argument to overcome this difficulty 

To start, consider for s E C+ the characteristic exponent of the rescaled number mea- 
sures 

ip(t,s)= / (1 - e~ sx )h(t,dx), <p*(s) = / (1 - e~ sx )n :¥ (x)dx 

J(0,oo) J(P,°°) 
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and let (fo(s) = <p(0,s). The Fourier-Laplace transform of the rescaled mass measures 
satisfy 



u(t,s) = d s ip(t,s) 

'(0,oo) 

with u (s) = u(0, s). By (1.12), 



c xn(t,dx), it*(s) 



(0,00) 



e sx xh*(x)dx 



(3.2) sup|F(t,x) -F*(x)\ < - sup 

x>0 x>0 



(u(t,a) — u*(cr))<icr 



+ 



24 

7tT' 



1. Evolution of characteristics: The evolution of ip and it is given by 
(3.3) d t <p + (2s - ip)d s cp = <p 



(3.4) 



dtu + (2s — (p)d s u = —u(l — u). 



As with the constant kernel, these equations can be derived from (1.3), or from (1.1) if 
the initial number measure has a density [7]. Fixing an initial point so £ C+, the unique 
solution to (3.3-3.4) is given in terms of the characteristics s(t;0, sq) which satisfy 



(3.5) 



als 

— = 2s-(p, s(0;0,s ) = s G C+. 



Along these characteristic curves we have that 



dtp 
~dt 



tp, 



du 
~dt 



-u(l — u) 



so that by integrating: 

(3.6) <p(t, s) = eVo(so), u(t,s) = e~ t 



uo{s ) 



l-«o(ao)(l-e-*)' 



These solutions are analytic, as shown in [6]. The explicit form for the characteristics 
so H- s, which depend on the initial data through (po, is then 

(3.7) s(t; 0, s ) = e 2i (s - M*o)0- ~ e "*))- 

This can be combined with the expression (3.6) for ip to give 



(3.8) 



<p(t,s) 



s 1 - (<po(s )/s )(l - e *) 

Note that since |«o(so)| < 1 an d I (so) I < |so| horn (3.6) and (3.8) we have 
(3.9) \ip(t,s)\<\s\, \u(t,s)\<l 

for s(t; 0, sq) such that sq G C + and t > 0. 
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The Fourier-Laplace transform 
(3.10) u*(s) 



1 



of the self-similar solution can also be stated in terms of characteristics. In contrast to the 
to the constant kernel case, these characteristics differ from those for the solution u(t, s) 
with initial data no- Fixing Sq E C, this distinct set of characteristics Sq i— >■ s* satisfies 

(3.11) s*(t; 0, «S) = e 2t (s* - p,(*5)(l - e"*)) 



with (/3*(s) = \/l + 2s — 1. Since u^s) is a stationary solution along these characteristics, 
we necessarily have that 



(3.12) 



it*(s) = e" 



u*(«o) 



l-«,( S 5)(l-e-*)- 
Estimates near the origin: We use the additional finite moment assumption 



(0,oo) 



x 3 no(a!x) 



(0,oo) 



x 3 no(dx) < oo 



along with 



/ x 3 h*(x)dx = / x 3 n*(x)dx 

J(0,oo) J(0,oo) 



(0,oo) 



to approximate uq and ti* near the origin. By (1.17), 



M*)-(i-«)l<f M 2 , 



|«*(s) - (1 - < -|s| 



(3.13) 

Since no = d s <po implies 

<Po( s ) = s u (/3s)dp, 
Jo 

together with an analogous expression for 92* these estimates give that 



s G 



(3.14) 



^o(s) 



1 . 

S ~2 S ' 



~ 6 |S| ' 



1 - 

s s' 

2 



1, 

< 5 M 



s e C_ 



This yields a bound on the distance between the characteristics s, s* starting at ctq G C+: 
\s(t; 0, cj ) - s*(t; 0, <r )| = e 2 Vo(^o) - <^*(<t )|(1 - e - *) 

(3.15) 



<(J + ^)(l-e-VVo| 3 



With (3.9) we also have that 



(3.16) 



s*(t;0,cJo)| = e M \a - <p*(a )(l - e"*) < 2e 2 > |- 



Finally, we estimate u(t, s) itself near the origin using the time-dependent moment 



m 3 (t) 



x s h(t, dx) 



(0,00) 



(0,00) 



x 3 n(i, dx). 
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It is straightforward from (1.3) (or (1.1) if n(t, dx) has a density) that dtrh^{t) = 3 — rh^[t), 
so 

m 3 (t) = 3 + (/i 3 -3)e-*. 
Then (1.17) implies that for all t > 0, 

(3.17) |u(t,s) - (1 - 8)| < — ^— \s\ < —\s\ . 

3. Propagation of estimates: Upon inspection of (3.7) and (3.11) we see that the natural 
growth rate of both sets of characteristics is e* near the origin, and e 2t at an 0(1) distance 
away from the origin. Explicitly, for small 5 > note that 



s(t; 0, iS) = e 2t (i5 - <po(iS)(l - e _t )) « e 2t (iSe' 1 + ^(i5) 2 (l 



Bending this point onto the imaginary axis, define 

(3.18) T{t) = e 2t Ue' 1 + U 2 (l - e"*) 

We will bound the r.h.s. of (3.2) with this value for T, where 5 = 5(fj,<2,) > is a length 
scale characterizing the region about the origin where the difference between no and u* is 
small by (3.13). To begin, let p 3 = max{// 3 ,3} and define 



(3-19) = A/I + K7 + ? 



2 V4 6 

so that (| + 5 = |(2 + (5)" 1 . Note that if /i 3 -)• 00 then 5^0. 

^. Backward characteristics: The main difficulty in the proof arises from the fact that 
characteristics are no longer rays as in the constant kernel case. Let f2j an d £1% be the 
images of C + under the maps sq h-> s(i;0, sq) and Sq i-> s*(i;0, s ). In [7], Lemma 3.3, it 
has been shown that the images I\ = and T* = c?0* of the imaginary axis under these 
maps are curves which pass through the origin and otherwise lie in the closed left half plane 
<C_. This is easily seen by taking the real part of (3.7) with s = ik and noting that 

Re (fo(ik) — I (1 — cos(kx))no(dx) > 0, k 7^ 0. 

J(0,oo) 

Furthermore, it has been shown in [7] that T t and Tj can be considered graphs over the 
imaginary axis. We will use these properties later in the proof. 

Define the backward characteristics so(t,s) and Sq(£, s) as inverses of the forward map- 
pings so ^ s(t;0,so) and Sq h-> s*(t;0, Sg), respectively. Explicitly, for and G O*, 

(3.20) s (t,w) = s(0;t,w), s (i, w*) = s*(0; t, w*). 

Also note that by (3.7) and (3.11), the Jacobian of the transformation sq 1— >■ s and s h> s* 
is 

(3.21) |i = e 2 *(l - n ( So )(l - e" 4 )), ^ = e 2i (l - u.(*S)(l - e~')). 
aso "Sq 
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We estimate (3.2) by considering separately the integral near and away from the origin. 
In doing so, we use the following convention for the remainder of the proof. The independent 
variable a is used solely for functions on the imaginary axis in the s,s*-plane, while <jq is 
used analogously in the sq, Sg-plane. 

4- Bounds for < |<r| < 5e~ l : Applying (3.13) and (3.17), we get that 



0<H<<5e-* a 



-(u(t, a) - u*(o))do- 



< 2 



\o~\ 



\u(t, a) — u*{p)\d\a\ 



(3.22) 



1 , ^3 \ 



<l 5 + f I* 1 " 



5. Bounds for 5e~ l < \o~\ < T(t): Since the flow of characteristics no longer leaves the 
imaginary axis invariant, we use a contour deformation argument. 
5a. Contour deformation: Take 

A 1 = i5e~ 2t , V = id 

to be the endpoints of a piece of the imaginary axis in the sq,Sq- plane (see Figure 3.1). 
Let AD and A*D* be the images of A'T>' under the mappings sq \— > s(t;0, sq) and Sq >— >• 
s*(t; 0, Sq). As discussed earlier, AD C Tt and A*D* C lie in the left half of the complex 
plane. 

The part of (3.2) that remains to be estimated is an integral over the imaginary axis in 
the s,s*-plane, bounded away from the origin, with endpoints 

B = i5e~ t , C' = iT{t). 

To do so, split the integral into its individual components 



(3.23) 



BC ° 



(u(t,a) — u*(o~))do~ = / u(t,a)da— / u^(a)da. 



BC ° 



BC a 



Since the integrands are analytic away from the origin, we can deform the contours by 
applying Cauchy's theorem: 

f f e sS: f e s ^ f e s£ 

(3.24) / u(t,a)do~ = / — u(t,s)ds+ / — u(t, s)ds + / — u(t,s)ds, 



BC 



AD 8 



BA s 



DC s 



and 
(3.25) 



BC a 



-u*(a)dcr 



A*D* s 



—u*(s*)ds* + 



BA* f 



-u*(s*)ds* + 



D*C s 



—u*(s*)ds* 



The error terms are integrals over contours that connect BC to AD and A*D*, respectively. 
We take AB and A*B to be straight line contours. Additionally, let DC and D*C be such 
that their backward images 



(3.26) 



V'C = s (t,DC), V'C*' = s* (t,D*C,) 
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A A 




Figure 3.1: Endpoints of contours on imaginary axis are B = i5e l ,C = iT(t) and A' = 
i5e- 2t ,V = iS. On the left, A = s(t;0,A),A* = s*(t;0,A),D = s(t;0,V),D* = s*(t;0,V). 
On the right, B' = s (t,B),B*' = a*(t,B),C' = s (t,C),C*' = s*(t,C). 



are straight lines. We will estimate the integrals over these regions in (3.24-3.25) at the end 
of the section. 

Making the change of variables s i— > sq and s* \— > Sq and using (3.6), (3.12), and (3.21), 
we have: 



(3.27) 



AD s 
-I 



-u(t, s)ds 



A'V 



A*D* s * 



■21 



s(t; 0,(7o 



s*(t;0,a ) 



u*(a ) I da . 



Using the basic inequality |ai&i — 0262) < |^i — &2 1 1 «i | + — a2||&2| ; we can bound the 
magnitude of (3.27) by the sum Ai + A2, where 



(3.28) 



Ai 



s 

6e~ 2t 



e 2t e 2t 

W ((T ) -W*(<7o) 

S S* 



\e s x \d\a \ 



(3.29) 



A, 



JSe- 2t 



sx s x\ 

e — e 



e 2t 



d\a \ 



We now evaluate each of these terms. For Ai the argument is nearly identical to that 
presented in Section 2. 

5b. Bounds for A\: Since s* £ C_, \e s x \ < 1. Using the explicit form (3.7), (3.11) of 
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the characteristics, we rearrange terms in the integrand of (3.28) to get 



& 2t e 2t 

— u (ao) -w*(cro) 

s s* 



(3.30) 



(TO 



TO 



a h_ i£s^22l(l _ e -t) 



10(0-0) 



(1 - er*) 



Notice that here </>o(c?"o)/co and </?*(<7o)/<7o play the same role as uo(o"o) and n^(o"o) do in 
(2.17). Using that < 1 and the estimates (3.15), (3.16), the numerator in (3.30) satisfies 



e~ 2t s* 



CO 



«o(°o) 



CO 



U*(CT ) 



< 



e 2i s* e 2t s 



CTQ 



(70 



K(e>o)| + K(o"o) - «*(o-q)| 



e" 2 V 



(3.31) 



<5 + ?)M' 



To estimate the denominator of (3.30), we replace </?o(co) and </?*(o"o) by 2co/(2 + cr )- 
Since 



2s 



2 + s 

with (3.14) and (3.19) this implies that 



1 , 

s s' 

2 



1, 

< - s 



s 6 



i-^M(i- e -*) 

O"0 



> 



2e-* + a 



2 + a 



4 + T' 1(701 



(3.32) 



> 



|o-q| 
2 + 5 



4 + 6 



kol 2 > 



2 2 + 5 



The same argument holds with ipo replaced by ip* and /U3 replaced by 3: 



(3.33) 



1 



<£*(o"o) 



O"0 



(1 - e-«) 



> 



1 kol 



2 2 + 5 



Therefore, using (3.33-3.32) and (3.31) in (3.30) yields 

\e s * £ \d\a \ 



<5e- 2t <|cro|<<5 



e 2t e 2t 
^o(o"o) -U*((To) 

s s* 



< 16 (\ + f ) (2 + Sfe-t f ^-d\a \ 

V 2 6 / he- 1 * F0| 



(3.34) 



32 



5c. Bounds for A%: Now we consider A2. Using |u*| < 1 and (3.11) we can rewrite this 



as 



(3.35) 



A 2 < e 



-/ 



<5 |g s ^ e s *^l 



5e~ 



y*(p'o) 



(1 - e -*) 



-rf|cr |. 
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For the term in the denominator, we use the estimate (3.33). Then, 

1 



se- 



mi 



1 



< 



y*(o~o) 



(1-e- 



d\a \ 



TO 



1 

Se~ 2t l°0 



■d\ao\ 



(3.36) 
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5 + f ) 



5d. Bounds for Ai + A2: Putting together (3.34) and (3.36), we have a bound for 
(3.28-3.29): 



(3.37) 



( \ + y ) 



6. Error terms: Lastly, we bound the two error terms in (3.24) obtained from the 
contour deformation argument. The two remaining error terms from (3.25) are bounded by 
identical arguments. 

First we consider the term 



(3.38) 



[ — u(t,s)ds < \BA\sup ((si" 1 ) 
Jba s ba 



Here we have used that for s G Tt, both \e sx \ and |u(t, s)| are bounded by 1. From the 
explicit expression for the endpoints we have that 

\BA\ = \e 2t {i5e- 2t - ^ (^e~ 2 *)(l - e - *)) - We - *|. 

Substituting (po(i8e~ 2t ) with i5e~ 2t and using the bound 

(3.39) 



1 



l^o («o) - s \ < ^kol 2 , s G C + 



derived from (1.17), we find that \BA\ < \b 2 e~ 2t . 

For the remaining term, write s G BA as an interpolation of the endpoints: 

s = p(e 2t {i5e- 2t - ^{i8e- 2t ){\ - e"'))) + (1 - P)i5e~\ G [0, 1]. 

Again, substituting ipo(i5e~ 2t ) by i5e~ 2t and using (3.39) implies that for s G BA, 

\ s \ > 6e -t _ ^2 e -2t > l §e -t 

1 1 - 2 ~ 2 

Putting this together with the previous bound, we conclude that 



(3.40) 



BA s 



-u(t, s)ds 



< Se-\ 
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Now we estimate the second error term. Changing variables from s to sq, 

e s(t;0,s )x Uq ( So ) 



DC s 



-u(t, s)ds 



vc 



*0 



1 



(1-e- 



-ds 



(3.41) 



< e'^VC'l sup (Isol -1 ) • sup 
VC" vc 



so 



The endpoints of the line segment D'C are i5 and so(t,iT(t)). By (3.7) 



-2/ 



iT(t) = s (t, iT{t)) - if (s (t, »T(i)))(l - e"*). 

Taking t — > oo, we see that so(t,iT(t)) converges to ttq G C + such that 

1 



(3.42) 



7I"0 - ^o(7To)- 



This implies that ttq, and hence the line T>'C, is bounded away from the origin. To see this, 
use (3.39) and take absolute values in (3.42) to get |-7To| > 5. Therefore, sup^,/^ | so I 1 < 
5- 1 = : d{5). 

We also have an upper bound for | vro | as follows. Since the mass measure xho(dx) is not 
concentrated at x = 0, \uq\ is bounded away from 1 in an annulus about the origin. Then 
there exists an e = e(5) > such that for all sq € D'C', 



(3.43) 



^o(so) 


<-L 


so 





M/3s )|d/3< 1-e. 



Taking sq = ttq and substituting this in (3.42) we obtain 



|7T 



< 



^o(tto) 



7TQ 



< 



2e' 



This shows that \T>'C'\ is bounded above by some constant C2(<5). In addition, (3.43) shows 
that 

^o(so) 



sup 

vc 



1 



so 



-(1-e- 



< £ 



Finally, putting these estimates together in (3.41) yields 



(3.44) 



DC s 



-u(t, s)ds 



< C 3 (<S)e-*. 



where C3 := CiC^e -1 depends on the initial data uq only through 

7. Rate of convergence: Using (3.22) and (3.37) in (3.2) with the bounds (3.40) and 
3.44) for the error terms, we obtain the desired result (3.1). 
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3.2 Monodisperse initial data with K = x + y 

As shown in [2], the explicit solution with K = x + y and initial data no(dx) = 5±(dx) can 
be given in terms of the Borel- Tanner distribution 

B(X,j) = X^ lJ -— ] , 0<A<1 

by 

oo 

(3.45) n(t,dx) = ^e _ *S(l - e'* , j)5j(dx), te[0,oo). 

3=1 

In the probabilistic setting, B(X,j) is the distribution of the total population size X\ of a 
Galton- Watson branching process with one progenitor and Poisson(A) offspring distribution. 
Note that 

oo 

E[X x ] = Y,jB(X,j) = (l-Xr 1 . 

3=1 

In similarity variables (1.9) the distribution function for xh(t,dx) corresponding to (3.45) 
is 

(3.46) F{t,x) = e~ 2t ^ e*(l - e^' -1 - 6 



l<j <e 2t x 



and the limiting distribution is 



(3.47) F.(x) = [ -}=r 1/2 e-y /2 d y . 



JO 

We consider now sup £>0 \F(t, x) — F*(x)\. Substituting Stirling's formula for the factorial 
in (3.46) and multiplying by an extra factor of 1 — e - *, define 

$(t,x) = e- 2 * Yl (1-e-^ ^ . . 

\<j<e*i V2ivf + 2e-i 

In addition, define the Riemann sum corresponding to (3.47) by 



l<j<e 2t x 

Each of the differences 





D 1 ~- 


= F(t,x)- 


$(t,x) 


(3.48) 


D 2 = 


= $(t,x) - 


$*(t,x 




D 3 ~- 


= $*(t,x) - 


-F*(x) 
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is estimated as follows. For D\, we can r eplace F(t,x) by (1 - e"*)F(t, x) since the dif- 
ference between these terms is 0(e ). Utilizing the following error estimate for Stirling's 
approximation [1] 

j\ = V2^j j+ h~ j e £(j) , < e(J) < — 



12j' 



we find that 



l<j<e 2t x ' 

= e~*(l - e~*)P(Xi_ e -t < e 2 *£) + 0(e~*). 

To estimate -D3, we use that the Riemann sum corresponding to an integral of a differentiate 
function leads to an error proportional to the grid size. Since F'J^x) G C 1 we have that 
D3 = 0(e~ 2t ). Finally, for the remaining term D 2 note that 

D 2 = e~ 2t Yl -^{e- 2t jr 1/2 e-^ 2t til 2 (e a ^ - 1 



/2tt 

l<j<e 2t x 

where 

«(*, J) = - log(l - e"*)) + e- 2 *j72 = -± e - 3 ' j + 0{e~ u ). 

Therefore, 

D 2 = V< [ X -}==yV 2 e-V 2 dy + 0(e- 2t ). 
Putting together these estimates for (3.48) we conclude that 

Bvp\F(t,x)-F*(x)\ = 0(e _t ). 

4 Rate of convergence for the multiplicative kernel K = xy 

Lastly, we consider the case of the multiplicative kernel: 

Theorem 4.1. Let no be a positive measure such that Jj , x 2 no(<ix) = Jj Q ^ x 3 ?io(dx) = 
1, and U4 = J^ ^ x 4 no(dx) < 00. With r = log(l — and n(t, dx) i/ie rescaled solution 
to Smoluchowski's equation with K = xy and initial data n (dx), let 

F(t,x)= / y 2 n(T,dy), F*(x) = / y 2 n^ 2 (y)dy. 
J(o,x] J{o,&] 

Then for r S [0, 00) , 

(4.1) sup\F(r,x)-F*(x)\ < CM(l + r)e- r 

where C{\i/±) is a constant that depends only on 
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4.1 Proof for the multiplicative kernel 

The rate for the multiplicative kernel can be recovered from that of the additive kernel 
by a classical change of variable due to Drake [4] (and discussed in [7] ) . With initial data 
no(dx), the measure-valued solutions n mnl (t, dx) and n ad(i (t,dx) to (1.3) with K = xy and 
K = x + y, respectively, are related by 

xn mul (t, dx) = (1 - t^n^irit), dx), t G (0, 1). 

Using the similarity variables 

x x 

x l = -T-77T) x 2 



xi(ty z A 2 (t) 

given by (1.7) and the rescaled number distributions (1.9), we have for r € [0, oo) that 

x 2 2 h mul (T,dx 2 ) = xin add (r,dxi). 

Therefore, 

F 2 (r,x)= f 2 h™\ T ,dy 2 )= y 1 n add (r,dy 1 ) = F 1 (r,x). 
Jo Jo 

Since (1.5) implies that the limiting distributions satisfy F*^{x) = F*,i(x), we combine this 
with the previous equation to get (4.1): 

sup|F 2 (r,£) - F, j3 (x)| = sup|Fi(r,x) - F %x (x)\ < C(/i 4 )(l + r)e~ T . 

x>0 x>0 

The constant C is the same as that obtained in the case K = x + y, with 113 replaced by 
/i4. Lastly, the convergence rate is nearly optimal for monodisperse initial data by making 
the change of variables to the additive kernel case. 
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